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Abstract 



We reconsider the valuation of barrier options by means of binomial trees from 
a "forward looking" prospective rather than the more conventional "backward 
induction" one used by standard approaches. This reformulation allows us to 
write closed-form expressions for the value of European and American put 
barrier-options on a non-dividend-paying stock. 
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I. INTRODUCTION 



Options are financial contracts which give the holder the right to buy (call options) or 
sell (put options) commodities or securities for a predetermined exercise (or strike) price by 
a certain expiration date Conventional European (American) options can be exercised 
only on (at any time up to) the expiration date. Since the option confers on its holder a 
right with no obligation, it should carry a price at the time of contract. It is the classic work 
of Black, Scholes and Merton which suggested a strategy for determining a fair price 
for the option in a risk-free environment. 

Closed-form valuation within the Black-Scholes-Merton equilibrium pricing theory [|2]|| 
is only possible for a small subset of financial derivatives. In the majority of cases one must 
appeal to numerical techniques such as Monte Carlo simulations, or finite difference methods 
and much of the effort in the field has been in developing efficient algorithms for numerically 
solving the Black-Scholes equation An alternative direction has been the evaluation 
of discrete-time, discrete-state stochastic models of the market on binomial and trinomial 
trees [f||| . Not only is this discrete approach intuitive and easily accessible to a less mathe- 
matically sophisticated audience; but it also seems to us to be a more accurate description 
of market dynamics and better suited for evaluating more involved financial instruments. 
Moreover, the few exact Black-Scholes results available can be recovered in the appropriate 
continuous-time trading limit. The main difficulty in pricing with binomial trees has been 
the non-monotonic numerical convergence and the dramatic increase in computational effort 
with increasing number of time steps |||7|]. For example, the state of the art calculations 
involve memory storage scaling linearly (quadratically) with the number of time steps, N, 
for European (American) options, while the computation time increases like N 2 in both 
cases [|J. 

In this paper we reconsider valuation on binomial trees from what we call a "forward 
looking" prospective: we imagine acting as well-educated consumers who attempt to elimi- 
nate risk and estimate the future expected value of an option according to some reasonable 
dynamical model. We will regard the movement of the price on the tree as a random walk 
(with statistical properties consistent with a risk-neutral world) with "walls" imposed by 
the nature of the option, such as the possibility of early exercise (American options) or the 
presence of barriers. The resulting mathematical formulation then has two conceptually dis- 
tinct components: the first ingredient is an explicit description of the possible "walls" . For 
example, in the case of barrier American options both the barrier and the "early exercise" 
surface need to be specified. The second step will be to compute the probability that the 
price reaches particular values at every accessible point on the tree. This involves counting 
the number of paths reaching that point in the presence of "walls", a somewhat involved 
but exactly solvable combinatorics problem. Once these two steps (specifying the walls and 
computing the probabilities) are accomplished the value of both European and American 
options, with and without barriers, can be written down explicitly. In an attempt to be 
pedagogical, we will limit ourselves to the simplest put options: European, simple American 
and European with a straight "up-and-out" barrier. Although the calculation can be simply 
extended to the barrier American option that discussion merits a separate publication. 

As far as we know, in the case of trees explicit formulas like the ones we are proposing 
exist in the literature only in the simplest case of conventional European options For 
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the more complicated case of American options, the main issues are best summarized in the 



last chapter of Neil Chriss' book ||10|| : "The true difficulty in pricing American options is 
determining exactly what the early exercise boundary looks like. If we could know this a 
priori for any option (e.g., by some sort of formula), we could produce pricing formulas for 
American options." Below we propose a solution to this problem in the context of binomial 
trees. Our formulation complements the earlier studies of American options in the limit 
of continuous-time trading [TT|^r5| which also focus on the presence of an early exercise 
boundary for the valuation of path-dependent instruments. The study of the continuum 
limit of our formulas is instructive and will be left for a future publication. 



II. BINOMIAL TREES 

To establish notation we begin by dividing the life of an option, T, into N time intervals 
of equal length, r = T/N. We assume that at each discrete time U — ir (i — 0, 1, 2, N) 
the stock price moves from its initial value, So = S(to = 0), to one of two new values: either 
up to Sou [u > 1) or down to Sod (d < 1) ||. This process defines a tree with nodes labeled 
by a two dimensional vector, (i,j) (i = 0, 1, 2, N; j = 0, and characterized by a 

stock price S(i,j) = SoU^d 1- ^ , the price reached at time ti = ir after j up and i — j down 
movements, starting from the original price So. The probability of an up (down) movement 
will be denoted by p u (pd = 1 — Pu)', and thus each point on the tree is also characterized by 
the probability, p{(l —p u )' l ~^ ) which represents the probability associated with a single path 
of i time steps, j (i—j) of which involve an increase (decrease) in the stock price. Computing 
the probability of connecting the origin with point (i, j) requires, in addition to the single 
path probability, a factor counting the number of such possible paths in the presence of a 
barrier and/or the possibility of early exercise. The calculation of this degeneracy factor 
involves the details of each financial derivative and it will be discussed in turn for each of 
our examples. 

The binomial tree model introduces three free parameters, u, d and p u . Two of these are 
usually fixed by requiring that the important statistical properties of the random process 
defined above, such as the mean and variance, coincide with those of the continuum Black- 
Scholes-Merton theory JT[|. In particular, 

p u u + (1 -p u )d = e TT (1) 
e rT (u +d) - ud -e 2rT = a 2 r, (2) 

where r is the risk-free interest rate, and the volatility, a, is a measure of the variance of 
the stock price. We are left with one free parameter which can be chosen to simplify the 
theoretical analysis; one might choose, for example, u = 1/d which simplifies the tree 
geometry by arranging that an up motion followed by a down motion leads to no change in 
the stock price. This condition together with ([!]) and (0) implies: 

u = e aVT (3) 

d = e-°^ (4) 
e rT — d 

Pu = -T- (5) 

u — d 
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We stress that Equations (1-5) are to be regarded as short-time approximations where terms 
higher order in r were ignored. 

With these definitions out of the way we can begin discussing the valuation of put options 
with strike price X and expiration time T. 



A. European Put Options 

The simple European put option is a good illustration of our "forward looking" approach. 
We are interested in all those paths on the tree which, at expiration time i = N, reach a 
price, S(N,j) = Sou^d N ~^ < X, for which the option should be exercised. That implies that 
j < j* = Int [ln(X/ Sod )/ \n(u/d)], where Int refers to the integer part of the quantity in 
square brackets. The mean value of the option at expiration can then be written as a sum 
over all values of j < j* of the payoff at j, X — SoU^d 1 ^ - ^, multiplied by the probability of 
realizing the price S(N,j) = Sou^d N ~^ after N time steps, P[N,j]. As already mentioned 
above, P[N,j] = N#[iV, j]pi(l — Pu) N ~^ ■, where KcfiV, j] counts the number of paths starting 
at the origin and reaching the price S(N,j) in N time steps. For the case of conventional 
European options this is just the number of paths of N time steps, with j up and N — j 
down movements of the price, and is thus given by the binomial coefficient, 

MJU= (?) («) 

The resulting expression for the mean value of the option at maturity is then discounted 
to the time of contract by the risk-free interest rate factor, e~ rT , to determine the current 
expected value of the option: 

Ve = e~ rT t[ N 3 ) riU ~ Pu) N ~> (X - S^d N ^) . (7) 

This expression is not new: it was first discussed by Cox and Rubinstein who also 
showed that in the appropriate continuous trading-time limit (r — > 0) (^) reduces to the 
Black-Scholes result [01. 



B. European Put Barrier Options 

We are now ready to extend (|7p into an exact formula for the mean value of an European 
put option with a barrier. Although our approach can be used for other barrier instruments, 
we consider the simplest case of an "up-and-out" put option which ceases to exist when some 
barrier price, H > So, higher than the current stock is reached. With the choice u = 1/d an 
explicit equation for the nodes of the tree which constitute the barrier can be written down: 

S(j B + l + 2h,j B + l + h) =S u jB+1+h d h , h = 0,l,...,h B (8) 

Here, js = Int [In (H/ S ) / ln(w)] defines the first point just above the barrier, (j B +l, j B +l), 
and Kb labels the last relevant point on the barrier corresponding to the time closest to the 
maturity of the option, i.e., h B = Int [(N — j B — l)/2]. 
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Since the probability that any allowed path starting with the present stock price, So, 
reaches an exercise price at maturity, S(N,j) < X, is still p{(l — p u ) N ^^ (with j < j*) the 
average value of the European barrier option can be written in a form similar to (]?]): 

V EB = e~ rT £ X EB [N,j] pi(l - p u f-3 (X - S uid N ~ j ) , (9) 

3=0 

where ^EB[N,j] is the number of paths N time-steps long involving j up and N — j down 
movements of the price excluding those paths reaching any of the points on or above the 
barrier (§). As we will explain below, ^E B [N,j] is given by 

Keb[N,j\ = ( N j) ~ ^^B\3B + l + 2hJ B + l + h]^.Zf B Z\Zf)^ ( 10 ) 

where the second term on the right-hand side represents the contribution from the unwanted 
paths which hit the barrier (^) before reaching an exercise point (N,j). 

To understand the form of the excluded contribution in (|H]) we first note that reaching 
the excluded region requires that the path hits the barrier at least once. One might think 
that the number of unwanted paths can then be calculated by (i) counting the number of 
paths connecting the origin to a given point on the barrier; (ii) multiplying this by the 
number of paths connecting that point on the barrier with the exercise point (N,j) ]rj| 
this includes all paths which wander into the above-barrier region; and finally (iii) summing 
over all points of the barrier (|Sp. However, a particular path reaching a given point on the 
barrier might have already hit any of the previous barrier points, and thus it would also be 
counted in the contribution in (ii) from all paths starting at the first barrier point reached by 
the particular path under consideration. Thus, summing indiscriminately over barrier points 
would lead to overcounting unless, in (i), we only include those paths which hit the barrier 
for the first time. In other words, (i) must only include paths starting from the origin which 
reach the particular point on the barrier without having previously visited any other barrier 
point. The number of such restricted paths (reaching the point (j B + 1 + 2h,j B + 1 + h)) is 
what we denoted by ^E B \j B + 1 + 2h,jB + l + h] in (pTOj) . Also note that the final sum over 
the length of the barrier is restricted to h < hu = min(hB,j — j B — 1) with j > j B + 1, 
corresponding to the fact that, in general, the exercise point (N,j) cannot be reached from 
all points on the barrier. This completes our explanation of (PIT). 

We are then left with computing ^ r £ B . From its very definition it is not hard to see that 
N^b[/i] = ^ r E S [jB + 1 + 2h,j B + 1 + h] satisfies the following recursion relation: 

(11) 

E^w^rj ). ^ ( i2 ) 

with the sum in (0) removing contributions from previously visited barrier points. Obvi- 
ously ^^[0] = 1 as there is a single path involving j B + 1 up moves connecting the origin 
with the point (j B + l,j B + 1) on the tree. 

To solve Equations (JTTJ) and flT2| ) we first combine the sum on the right-hand side of 
(|12|) with the term on the left and rewrite the resulting equation in the form of a discrete 
convolution: 



W£ B [h] 



j B + 1 + 2h 
jB + l + h 
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1=0 



2(h-l) 
h-l 



3b + 1 + 2h 



(13) 



where the boundary condition, ^£#[0] = 1, is already included as the h = contribution 
to (|TB|). Note that ( |TB| ) can be solved by standard Laplace transform (or Z-transform) 
techniques |L5|]. Since in applying these ideas to the more complicated American options 
we will lose the convolution form - the kernel will depend on h and I separately and not 
only through the difference, h — I - we will proceed in a more general way and stay in 
"configuration space" until the very end. 

We prefer to regard (|T3|) as a matrix equation of the form: 



t Tyres 



D 



EB- 



(14) 



Here II^b and Deb are hj^ + 1 dimensional vectors, with components Yl E B,h = ^ebV 1 ] 



and T>eb„ 



3b + 1 + 2h 
jB + l + h 



h = 0, 1, 2, /ijv/, and the {h M + 1) x (Km + 1) dimensional 



matrix, L E b, can be written as, 



j eb\ K i 



2{h-l) 
h-l 



9 (h-l). 



(15) 



Note that in (|T^) we have explicitly added a 9 function (9 (x) = 1 for x > and vanishes 
otherwise) to stress that L*eb is a lower triangular matrix with unity along and zeros above 
the diagonal. This simple observation allows us to rewrite (|T4|) in the convenient form, 



J EB 



h-l J 
solution for ([14]) 



{h M +l)x(h M +l) + Qeb- (16) 
for y > h M + 1; and [Qeb]^ = 
9{h — I — 1). The nilpotent property of Qeb allows us to write down the explicit 



where Qeb is a nilpotent matrix of order Q^^^ 
(2{h-l)\ 



U EB = [1 + Qeb] 1 D 



EB 



h M 

£( 

R=0 



•1) r QpdD 



EB EBi 



which in turn leads to the following formula for the value of the option, 
V EB = V E 

'N-j B 



V-res 
EB 



-rT 



V-res 
EB 

j* h M (j) 
j=j B +l h,l,R=0 \ J JB 

x p> u (l - Pu f-i (X - S u^d N -^ 



2h 
h 



Qeb 



h.l 



JB + 1 + 21 
JB + I + I 



(17) 

(18) 
(19) 



The lower limit, j = js + 1, on the external sum in (^) excludes all paths unaffected by 
the presence of the barrier; also we have explicitly indicated the iV and/or j dependence of 
the various quantities involved; and have separated out the contribution to K EB [iV, j] from 
unrestricted paths (the first term on the right-hand side of (|10D) which simply leads to the 
value of the European put option given in (^). 
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We expect that, since we have an analytical formula, we should be able to recover the 



exact solution of the continuum Black-Scholes theory for this simplest of barrier options [16 
as was already done for conventional European puts ||. Figure [I] shows the numerical 
convergence of the binomial value of a representative "up-and-out" European put option to 



its analytic value [16 



The same general idea used in the case of European barrier options will now be used to 
write down an exact formula for the price of a simple American option, regarding the latter 
as an option with an early-exercise barrier. 



C. Conventional American Put Options 

Using this view to valuate American options requires the knowledge of those points on the 
tree where it first becomes profitable to exercise the option. This set of points, parametrized 
as {hjx[i])i constitute the "early exercise barrier" (EXB). Determining the explicit form of 
the surface, j x [i], seems very difficult (if at all possible) as it already implies a knowledge 
of the mean value of the option at some finite number of points on the tree. In this section 
we show that there is a self-consistent exact formulation of the problem which proceeds in 
the following three steps: (i) we assume that the early exercise surface, j x [i], is given and 
compute an explicit formula for the value of the option at each point on the tree, f(i, j;j x [i}), 
which depends parametrically on j x [i]', (ii) the fact that early exercise at only occurs 
when X — Souid 1 '! > f(i,j',jx[i]) gives us an explicit formula for the EXB which corresponds 
to the strict equality, 

(■V " : ' ' ) = f(i,j x m .lr'. In, {./,:. /,./}. (20) 



[Note that, on the right-hand side of fl20"D we have not used f(i,j[i]) which might appear 
at first sight as a more natural choice for defining the EXB. As will become clear below, 
(p0|) is the simplest and most natural choice which resolves the ambiguity of defining f(i,j) 



away from points on the tree.] Finally, (iii) substituting the solution (EDI) into the formally 



exact valuation expression gives us the value of the option. Although this strategy leads to 
an exact solution of the price of an American option, explicit numbers require rather heavy 
numerical computations except in the simplest example of a straight EXB. 

Let us proceed in carrying out the program outlined above by assuming that the EXB, 
i.e., j x [i], is explicitly given. To begin our calculation we will need some very general 
properties of the barrier. These follow from two simple characteristics of early exercise: (i) 
if the point is an early exercise point, then so are all points "deeper in-the-money", 
f = 0; I; 3 ~ 1; an d (ii) if two adjacent points at the same time step, (i + 1, j + 1) 
and (i + 1, j), are both early exercise points so is the point (The latter property 

follows from a conventional "backwardation" argument |1J which indicates that the average 
expected payoff at discounted at the risk-free interest rate, is smaller than the actual 

payoff, thus making itself an early exercise point.) It is not hard to see that (i) and (ii) 
guarantee that the inner part of the early exercise region cannot be reached without crossing 
the EXB. Thus, if we define %a to be the first time for which early exercise becomes possible 
and parametrize the points on the EXB as (i = iA + h, j x [iA + h}) with h = 0,1,2, N — %a, 
it then follows that j x [m] = 0. Moreover, the structure of the tree ensures that j x [i] is a non 
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decreasing function of i; more precisely, for each time step, j x [i] either increases by one or 
remains the same. 

The formal expression for the price of an American option can be written down once 
one recognizes that once a path hits the EXB the option expires and thus any point on 
the barrier can be reached at most once. As a result, the value of the option is a sum of 
(appropriately discounted) payoffs along the barrier, weighted by the probability of reaching 
each point on the barrier without having visited the barrier at previous times. We can then 
write the expected value of an American option as: 

N-i A 

V A = e- r{iA+h)r ^ s [h}pt [iA+h] (l - Pu yA+h-MiA+h] ^ x _ SoU j x [iA+h] d i A +h-MiA+h}\ ; 



h=0 



(21) 



where N A es denotes the number of paths reaching the EXB in %a + h time steps without 
having previously visited any points on the barrier. 

The counting problem can be solved along similar lines to those followed in the case of 
European options: ^^ es [/i] satisfies an equation analogous to (0), namely, 

^[0] = 1 (22) 

*™ = ( ) " § *™ ( Ui, + % - 'L + 1] ) • h> - (23) 

where the first term on the right-hand side counts the total number of unrestricted paths 
from the origin to the point (i^ + h, j x [iA + h}) on the barrier, while the second term excludes 
those paths which, before reaching {%a + h,j x [iA + h}) visited any of the previous barrier 
points, (i A + l,jx[iA + l}), I = 0, 1, 2, h - 1 0. 



As in the case of the European barrier option ( |23| ) is rewritten as a matrix equation: 

L A Il A es = B A . (24) 

Here Ua and Y)a are N — i A + 1 dimensional vectors, with components n^ /, = W A es [h] and 

Da h = ( ■ \ A \i ] , h = 0, 1, 2, N — iA, and the (iV — %a + 1) X (N — %a + 1) dimensional 

\3x[i'A + n\ ) 

matrix, L^, takes the form, 

^,,= (^4:^ + ,,). «< k = 0,l,2,..., W -u. (25) 



Note that, in contrast to ( p~5| ) and (|T^), depends on the indices h and I separately; also, 



we have used the identities j x [i A ] = and j J = lj to incorporate the boundary 

condition, N^ es [0] = 1, in (^4|) in a symmetric way. As in (16), we can decompose L A as, 



M w = 6 hjl + [Q A ] h>l (26) 

h-l 

j x [i A + h] - j x {i A + I] 



[QaIm = ( , \ u^n)Hh-l-l), (27) 



S 



where Qa has nonzero elements starting just below the diagonal and it is thus a nilpotent 
matrix of degree N — i A (i.e., Q A ~ lA+l = 0). Thus, 



it 



N-i A 

[1 + Q^]- i D A = £ (— l) m Q™D j4 , 

m=0 

leading in turn to the final formula for the value of the option, 

3 — r(iA+h)T I i \m r/~\»ni I ^A I 



(28) 



V A = £ 



N—lA 

h,l,m=0 ' 



"(-i) m ra,„ 



(29) 



One last step is the determination of f(i,j), the value of the American put at every point 
on the tree which, in turn, will allow us to derive the equation for the EXB. This is 
easily done by simply translating the origin in (f29|): 



N-i A 

E e 

h,l,m=i—iA 



i A + l — i 
j x [iA + 



f(i,j)= E e-^ +ft -^(-l) m [Q: 

x pix[iA+h]-jQ _ p u yA+h-jx[iA+h]-i+j ^ _ g oU jxliA+h]^i A +h-j x [iA+h]^ 



(30) 



Together with (p0|) this then leads to the rather formidable-looking equation for the barrier 
height j[iA + k] at the (i^ + k)-th time step (k = 0, 1, ...,N — i^), as a functional of the 
barrier position at all future time steps before expiration: 



(X - S u ][iA+k] d iA+k ' ][iA+k] 



N-i A 

E • 

h,l,m=k 



-r(h—k)r i 



-l) m [Q 



/ - k 



A\h,l 



jx[iA + k] 



coefficient 



j x {i>A + I] ~ 3x\}A + k] J 

x pjx[tA+h]-j x [iA+k] ^ _ p^h-k-jx[iA+h]+jx[iA+k] 

x (X - S u j * [iA+h] d iA+h - jx[iA+h] ) 

= lnt{j[i A + k]}. (32) 

0|) j must be restricted to points on the tree as the binomial 
would be ill-defined - hence the choice (150).] Equations 



[It should now be clear that in 



j x [x A + 1} ~j[x A + 1} 

(|3T|) and (|32"| ) for the boundary together with the formula for the value of the option, (p9|). 
constitute an exact pricing strategy for a conventional American put. A similar formula for 
an American put with an "up-and-out" barrier will be discussed in a future publication. 

It is instructive to consider Equations (|29|) , (|31|) and (|32|) in the explicitly solvable case 
of a straight barrier. We begin with the observation that, at expiration, k = N — i A) ( plf ) 



reduces to the equation for j* = Int[ln(X/S , c/ Ar )/ \n(u/d)}, already defined in the case of the 
European option, and thus, the barrier goes through the point (N,j*). Moreover, starting 
from the exact point (N,j*) on the barrier and decreasing j x [i] by one with each backward 
time step we reach i A = N — j* along the straight line, j x [i] = i — N + j* . Recall that, 
since with each increasing time step, j x [i] either increases by one or remains the same, this 
straight line represents a lower bound for the early exercise barrier. 
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For this straight barrier ( p8|) and fl29|) reduce to, 

p s t ra ia ht = e -r(N-r)r ^ _ p jN-j> ^ _ g^N-j* 
h=l 



+ £ e~ r{N ^ +h)T ^f mght [h] p h u (1 - pvf-'* (X - S u h d N ~ jt ) (33) 



with 



*r wht = r i 1 ) ( " i 'r x i (34) 





We expect that the result for the true barrier should approach the straight line formula 
for coarse enough time steps, r > N —j* — ia, (where this %a is the first time of early exercise 
in the limit of continuous-time trading). 

III. CONCLUSION 

We have presented a scheme for pricing options with and without barriers on binomial 
trees. To the best of our knowledge ours is the first explicit derivation of exact formulas 
treating barriers on binomial trees. It is our expectation that in the limit of continuous- 
time trading we should be able to recover the few exact results available in the literature, 
especially for American options [[j3)[0l • We also hope that our explicit formulas may provide 
a framework for improving the efficiency of numerical computations. 
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FIGURES 




FIG. 1. Convergence to analytic value |16| of a three- month European "up-and-out" put 
option on a non-dividend-paying stock as the number of binomial intervals N increases. The solid 
line joining the squares is a guide to the eye. The stock price So is 60, the risk-free interest rate r 
is 10%, and the volatility a is 45%. The barrier level H is at 64. The analytic value (continuous 
line) is 2.524. As a point of reference, we give the European put option value (i.e., H —* oo) 4.6 
after the Black-Scholes solution (dashed line). 
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